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Abstract 

Cosmological dynamics of scalar field with a monomial potential (fP with a general 
background equation of state is revisited. It is known that if n is smaller than a 
critical value, the scalar field exhibits a coherent oscillation and if n is larger it obeys 
a scaling solution without oscillation. We study in detail the case where n is equal to 
the critical value, and hnd a peculiar scalar dynamics which is neither oscillating nor 
scaling solution, and we call it a pseudo scaling solution. We also discuss cosmological 
implications of a pseudo scaling scalar dynamics, such as the curvature perturbation 
and the domain wall problem. 



1 Introduction 


In the early universe, a scalar field often oscillates coherently after inflation [^. Coher¬ 
ent oscillation of a scalar field may have various cosmological effects such as generation of 
curvature perturbation, dark matter, baryon asymmetry and unwanted relics. 

In a simple model of a real scalar field 0 with a canonical kinetic term and a scalar 
potential of C oc 0", it is known that the scalar field stops at the initial position 0j until the 
Hubble parameter becomes equal to the effective mass of the scalar. After that, 0 begins a 
coherent oscillation for n < ric with ric being the critical exponent, which depends on the 


background equation of state (Eq. (2.10)). In this case, the mass scale of 0 becomes larger 


and larger compared with the Hubble scale and such a coherent scalar dissipates its energy 
due to the Hubble expansion. The amplitude of the coherent oscillation d) decreases as 


$ oc 


( 1 . 1 ) 


with a{t) being the scale factor of the universe. Such a scalar field can act as a curvaton [^, 
which generates the primordial density perturbation of the universe [^. In particular, cur¬ 
vaton models with n > 2 were studied in some literature [5]-[^. 

On the other hand, for n > ric, the scalar field never oscillates and it just follows a 
so-called scaling solution 

0 cx . (1.2) 

This is an attractor and it was pointed out that such a scalar field cannot generate curvature 
perturbation since the initial fluctuation of 0i is erased during the scaling regime |^. 

In this paper, we revisit the scalar dynamics for n = ric. In Ref. it is stated that in 
this critical case, the scalar field is not attracted to a scaling solution but it is something like 
between a scaling and oscillation. Nevertheless, they say that for n = ric both the oscillating 


solution (1.1) and scaling solution (1.2) show the same behavior with the cosmic expansion 


and hence the curvature perturbation also vanishes. However, we will see that the actual 
dynamics of the scalar field for n = ric is far from both the oscillating solution (1.1) and 


scaling solution (1.2). Even averaged features of the scalar dynamics cannot be described by 


these solutions. Instead, a peculiar feature of the scalar dynamics shows up, which we call a 
“pseudo scaling solution”. We will show it both in a numerical and analytical way and see 
that the curvature perturbation does not vanish in this case. 

The pseudo scaling solution can have a variety of physical implications. As we noted 
above, the curvature perturbation can be generated by a curvaton even if it follows the 
pseudo scaling solution, and hence it may offer the seed of the large scale structure of our 
universe. Another example is a symmetry breaking field such as the Peccei-Quinn field [^. 
In the supersymmetric (SUSY) axion model, it is possible that the Peccei-Quinn field follows 
the pseudo scaling solution in the early universe. As we will see later, we can naturally solve 
the isocurvature perturbation as well as the domain wall problem in such a case. The pseudo 


Note on terminology: in this paper we just call Eq. (1.2) scaling solution. We do not regard an oscillating 


solution as a scaling solution even if its amplitude satisfies a simple power law like (1.1). 
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scaling behavior also has some implications on the baryon nnmber asymmetry produced by 
the Affleck-Dine mechanism 10 . Therefore, it is worth studying detailed features of the 


pseudo scaling solution. 

The organization of this paper is as follows. In Sec. we discuss the dynamics of the 
scalar held with a monomial potential in a general background equation of state. Here 
the “pseudo scaling” scalar dynamics is extensively studied. Also, some comments on the 
Hubble mass term are given there. Cosmological implications of the pseudo scaling solution 
are explained in Sec. We discuss the curvature perturbation and the domain wall problem 
there. The last section is devoted to conclusions and discussion. 


2 Scalar dynamics with monomial potential 

2.1 Scaling and oscillating solution for n ^ ric 

Let us consider the Lagrangian 


£ = ( 2 . 1 ) 

where the scalar potential is given by the monomial one 

V{4>) = -r, ( 2 . 2 ) 

n 

where n > 2 is an integer]^ The equation of motion of 0 is given by 

0 + = 0, (2.3) 

with H being the Hubble parameter. We take H = p/t and assume 3p — 1 > 0. In the 
matter (radiation)-dominated (MD (RD)) universe, p = 2/3 (1/2). A more general equation 
of state may be obtained in the inhaton oscillation dominated universe (see App. 0 . Let us 
dehne 


(p = 0 , 

s = \n{t/ti), 


(2.4) 

(2.5) 


with a{t) being the scale factor of the universe and ti the initial time. We have taken 


a{ti) = 1. Then the equation of motion becomes 11,12 


//, /^q n + 2\ 2 , V n-i_n 

(p + ^3p — ^ _2 j ^ + p (p + — Oi 


( 2 . 6 ) 


For an odd n, 0 should be interpreted as |0|. 
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( 2 . 7 ) 


where the prime denotes derivative with respect to s and 


2 _ 2(6p — 3pn + n) 


{n - 2y 


The equation (|2.6|) represents the motion of the rescaled held p in the effective potential 

(2.8) 

In these variables, p has a minimum at (pmin where 


T/ ^ ^ ^ 2 2 , V 

^efr(y^) - ip + ^^2 


Thus if 


n > Tlr = 


min 


6p + 2 

3p- 1 


if 

-p 2 i 72 /^p 2 ^ 1 /(n- 2 ) if ^ 


> 0 
2 < 0 . 


(2.9) 


for p = 1 
for p = 2/3 
for p = 1/2 


■f-)- p = 


?7-c + 2 

3(ne-2)’ 


( 2 . 10 ) 


the sign of the friction term is positive, and hence the oscillation of p damps and it relaxes 
to the minimum of the potential prmn- 

Let us suppose the condition (2.10) is satished so that there is a nontrivial minimum of 
the effective potential and p relaxes there, ^ This represents the scaling solution; in terms 
of the original held 0, 

(j){t) cx oc (2.11) 

In this case, independently of the initial condition, 0 hnally follows the scaling solution 
without oscillation. Note also that the scaling solution implies that the ehective mass scale 
of the scalar decreases in proportion to the Hubble parameter: rrigg = cx H^. 

On the other hand, if the inequality (2.10) is inverted, the friction term of (2.6) has a 
negative sign, and hence p never relaxes to its potential minimum. This corresponds to an 
oscillating solution. This behavior is also understood as follows. Assuming an oscillating 
solution for 0 with potential ~ 0"', we obtain 

$(t) oc oc t-6p/(^+2)^ (2.12) 


by making use of the Virial theorem with $ being the oscillation amplitude of 0. Thus the 
effective mass of 0 scales as mgff cx <l)(”'“ 2)/2 f-3p(n-2)/(n+2)_ jf continue to exceed 

the Hubble parameter H ~ t~^, the oscillation survives. This condition is written as 


3p(n — 2) 
(n + 2) 




This is just an inversion of the condition (2.10). 


n < 


6p + 2 

3p- 1‘ 


( 2 . 13 ) 


Note that < 0 <t=> ric + 2 < 2n, and 
nontrivial minimum of the potential if n > Uc 
potential (see Eq. (2.35)). 


Uc > 2 for p > 1/3. Thus, < 0 is ensured and there is a 
. It does not always hold if there is a Hubble mass term in the 
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2.2 Pseudo scaling solution for n = ric 


We have seen that a scalar held approaches to a scaling solution for n > ric and oscillating 
solution for n < independently of the initial condition. However, the case oi n = ric is 
more nontrivial. In this case, the friction term of (2.6) vanishes: 






(n + 2)2 A 


(n-2)2^ 9{n-2yHf 




n—1 


= 0 . 


(2.14) 


In terms of (^, it does not universally approach to some hxed solutions and continues to 
oscillate without damping. The resulting dynamics signihcantly depends on the initial con¬ 
dition of 0. If the initial condition is chosen so that ip oscillates around (pmin without crossing 
(p = 0, 0 also does not cross 0 = 0. The resulting dynamics of 0 looks like neither scaling 
nor oscillating solution; we call this a “pseudo scaling solution”. Hereafter we consider the 
case oi n = ric- 

Before solving the equation of motion, we must specify the initial condition. The natural 
initial condition is set just after inhation as 0(ti) = 0* and 0(tj) = 0j ~ — 00 . 1 /( 00 /3iJj, 
where Hi and are the Hubble parameter and the cosmic time at the end of the inflation, 
respectively. We assume 0j > 0 without loss of generality. Here <pi should satisfy 


mlM) = ;= H?. 


(2.15) 


since otherwise 0 rolls down the potential during inflation. This condition roughly means 
ifi < ipmin- Therefore, in terms of ip, it is initially located close to the origin and rolls down 
the effective potential (2.8) toward the minimum. In particular, we consider the case 

iPi ’Pvam- (2.16) 


By noting that the initial velocity of ip is given by 




6 

n -I- 2 



we obtain the conserved energy density of ip\ 


Pip 


1 

2 


Pi T 




n + 2 

3(n - 2fHi 


^20j0j + 


n + 2 A 07 A 

3n Hi ) 


n + 2 X(j)^ 

9n{n — 2) Hf 


(2.17) 


(2.18) 


Here we have used the approximation |0j| |77i0i| which is valid as long as ipi Pmin- 

Since the energy density is negative, ip does not overshoot the origin (^ = 0, which means 
that 0 does not oscillate around 0 = 0. The rescaled held ip just oscillates around ip min 
without damping]^ 

This observation crucially depends on our assumption on the initial condition 3Hi^i ~ — Bh Although 
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Figure 1: Time evolution of 0(t) for n = ric for the potential (jf' for p = 2/3 (n = 6) in the 
left panel and p = l/2(n = 10) in the right panel. We have normalized 0 by the Planck 

mass Mp. The time is normalized by where is a natural mass 

scale of this system. The red solid line represents a numerical solution, which corresponds to 
a “pseudo scaling” solution, while the green dashed line shows a hypothetical scaling solution 
for representative purpose. 


The dynamics looks so simple in terms of the rescaled held p. However, this “pseudo 
scaling” solution is more nontrivial in terms of the original held 0. Let us see Fig.where 
we have numerically solved the equation of motion of 0 with initial conditions iLj = 1, 0* = 
O.2(pmin,0j = —V'{(j)i)/3Hi for p = 2/3 {n = 6) in the left panel and p = 1/2 (n = 10) in 

the right panel. We have normalized 0 by the Planck mass Mp. The time is normalized by 

^-2\V2 . 




\ I ‘2i 

, where ^ natural mass scale of this system. The red solid line 

corresponds to a numerical solution, which corresponds to a pseudo scaling solution, while 
the green dashed line shows a scaling solution 0 cx t“ 2 /(n- 2 ) representative purpose. For 
comparison, we also show time evolution of 0(f) for cases with n = 8 (scaling), n = 6 (pseudo 
scaling) and n = 4 (oscillating) for p = 2/3 (n^ = 6) in Fig. It is clear from this hgure 
that the pseudo scaling solution resembles neither the scaling solution nor the oscillating 
solution even qualitatively. To understand this behavior, let us solve the equation (2.14) 
more carefully. 

Our goal is to understand main features of the pseudo scaling solution and estimate how 
the periodicity depends on the initial condition 0j,^^ Here it is important to note that the 


it holds during inflation, the relation is modified after inflation depending on the effective equation of state 
of the inflaton oscillation dominated universe as Hcf) — —pF'/(3p+ 1) for \V" 1. The approximation 
~ —V as an initial condition is valid for inflation models in which the inflaton mass scale is much 
larger than the Hubble scale, such as new inflation. We have numerically checked that is indeed negative 
for such a case. On the other hand, for chaotic inflation, our approximation is not justified and Pip may be 
positive. However, the dynamics is almost identical to the case of negative p^ after 0 is replaced with |0| 
even if 0 once overshoots 0 = 0 as long as (j)i ^ OiV'/Hi) is satisfied. 


We can formally solve Eq. (2.14) (see App. fl^, but the following discussion is enough for our purpose. 
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Figure 2: Time evolution of 0(t) for n = 8, 6 and 4 for p = 2/3 [ric = 6). The green dashed 
line is the solution for n = 8 (scaling). The purple line corresponds to the solution for n = 6 
(pseudo-scaling), and the cyan line is the solution for n = 4 (oscillating). 


dynamics of cp is dominated by the region p p 
our interest 


#6 


since \p^\ -C |14((^min)| in the case of 


It means that we can almost neglect the second term of Eq.(2.8) in order to 


understand main features of the dynamics. In the region p p 
is approximated as 

„ 4 

^ - 1 -= 0 - 

[n — ly 

This has a general solution consisting of the growing and decaying mode, 


1 , the equation of motion 


(2.19) 


p{s) = A exp 


2s 

n — 2 


+ B exp 


-2s 


n — 2 


( 2 . 20 ) 


where A and B are determined by the initial conditions. In terms of 0, this solution corre¬ 
sponds to 

/ , N 4/(n-2) 

c^(t) = A + B(^Aj . ( 2 . 21 ) 

Let us divide the time into the first phase ti < t < ti, the second phase ti < t < t 2 and 
the third phase t 2 < t < t^. Here, we define ti as the time at which the motion of 0 starts. 
The time when the motion of (j) stops is denoted as t 2 . Some time after t = t 2 the held 0 
again starts to move, and we dehne that time as (see Fig.[^. The p^ term in the effective 
potential becomes important only for some short periods at around t = ti,t 2 and t^. 

This also holds for other initial conditions as long as the relation 0^ ^ (P {V'/Hi) is satisfied. 
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First phase In the time interval ti < t < ti, A and B are determined by the initial 
conditions: 


A 


1 

2 





( 2 . 22 ) 


This solution breaks down at s ~ Si (t ~ ti) where becomes close to (fimin- After that 
passes through (pmin and goes back to ~ (fimin within order one unit time in terms of s. We 
can estimate si and ti as 


n 


Si ~ 


In 


iPr 




O tl ~ ti 






(n-2)l2 




n—2 


1/2 


^eff (0i) 


(2.23) 


Thus the time ti is roughly the epoch at which the Hubble parameter becomes equal to the 
effective mass of 0. It is obvious that 0 remains almost constant during this period. 


Second phase Next, let us consider the interval ti < t < t 2 - Just what we have to do is 
to determine the constants A and B in (2.20) with appropriate initial conditions. We set the 
initial condition aX ip = ipj (and s = sj), with ipj taken somewhere between ipt and praim say, 
Pj ~ 0.1</5min- The precise position of pj is not important since the dynamics is dominated 
by the region p -C Pmin as we noted before. Since is conserved during the whole process, 
we can express p'^ as 


p. = 


2p^ + 


[n 


2) 




1/2 


where 


e = 


rA — A Xp'^ 
36n . 


n — 2 



(2.24) 

«1 

,2,^n-2 ^ 

(2.25) 


T 3 \ ' 'r^'rmin 

The sign of p' is flipped because p changes the direction of motion due to the higher order 
term in Kfr(/5) at t ~ ti. Then, dehning Aj and Bj as the constants A and B in Eq. (2.20) 
at s = Sj, we obtain 


A/ ^ ]^PjP 


Bj ~ Pj, 


(2.26) 


hence Aj <C Bj. Then the end of the second period, s = 52, is the time at which the decaying 
mode ~ iJ^.e- 2 (s 2 -si)/(ri.- 2 ) ]3eQomes equal to the growing mode 
obtain S 2 (^ 2 ) as 


p^.^2[s2-sp/(n-2)^ 


n 


S2 


C . 

— 


-2 2 


O t2 — Cti 


Pr 


Pi 


n(n—2)/4 


(2.27) 


where in the last equality we have approximated tj ~ ti and pj ~ Praia- Here we have 
included an 0(1-10) constant c which comes from the uncertainty of pj as well as effects of 
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the 99 " term. We have numerically checked that the dependence of ^2 is precisely given 
by Eq. (2.27) for (p* 99min, or mes{(l>i) Hi. It should be noticed that ^2 is orders of 
magnitude longer than the inverse mass scale of the scalar held. From (2.21), we can see 
that (f) scales as oc (Juring this time interval, so we have 


0(t) ~ ipi 


/ X 4/(n-2) 

(f) 


for ti < t < t 2 , 


and eventually at t = ^2 h becomes 


0 (^ 2 ) ~ <Pi 



(2.28) 


(2.29) 


Note that during this regime, (j) satishes (j) + 3H(j) ~ 0, so the potential term becomes 
irrelevant and the effective mass scale mes{(p) becomes smaller and smaller than the Hubble 
parameter. The energy density of 0 is dominated by the kinetic term and scales as oc 
during this regime. 


Third phase Lastly, let us consider the third period t2<t < fg. This is just a continuation 
of the second period, where the growing mode dominates over the decaying mode. Thus the 
solution is 

2 


This is valid until it becomes close to (prain- 


n 


2 


(2.30) 


estimated as^’^ 


n 


S3 - S2 ~ 


2,2 / (^min 

-— In -, o fs ~ ct 2 

4 e 


Thus the end of the third period, S 3 (fs) is 

n(n—2)/4 


(2.31) 


Pi ) meff(0(t2))’ 

Therefore it is seen that the relative time interval ^ 3/^2 is equal to ^ 2 /^- During this time 
interval, 0(t) = 0 (^ 2 ) remains almost constant as is clear from (2.21). 


After f = fs, we have exactly the same dynamics as the second phase and the following 
third phase since there is no damping in the equation of motion of p. This one “pseudo¬ 
oscillation” consisting of one set of the second and third period takes much longer time than 
the inverse mass scale of the scalar held and even than the Hubble time scale. In terms of 
0, the held value reduces by an amount of e ~ {pilPmix^'^ during this one pseudo oscillation. 
The scaling of 0 consists of 0 = const regime and 0 oc ^“^(>^- 2 ) regime, both are completely 
diherent from the scaling solution oc t“ 2 /(n- 2 )_ scaling solution is obtained only after 
averaging the pseudo oscillation, but the period of the pseudo oscillation is much longer than 
the Hubble time scale, hence using the scaling solution would lead to qualitatively wrong 
conclusions on the physical implications such as the curvature perturbation. 


The constant c here is exactly the same as that in Eq. (2.27) since the dynamics of y is periodic. 













Finally, we make some comments on qualitative differences between the scaling solution 
for n > Uc and the pseudo scaling solution for n = ric- In the case of the scaling solution, the 
kinetic term, the Hubble friction term and the potential term in the equation of motion are 
all comparable. Therefore, if the potential is not monomial, the motion of 0 soon deviates 
from the scaling solution once other terms than the original one dominate the potential. On 
the other hand, in the case of the pseudo scaling solution, the kinetic term and the Hubble 
friction term dominate over the potential term. This can be seen by the fact that (j) = const 
and (f) oc the solutions of the equation 


(j) + 3H(j) = 0, (2.32) 

for p = (n + 2)/3(n — 2),^® Thus, even after the potential is dominated by other terms than 
the one with n = ric, the motion of 0 still follows the pseudo scaling solution for some period 
until the potential term becomes comparable to the kinetic and friction termsp^lWe will see 
this in more detail in Sec. |3l 


2.3 Comment on the case of Hubble mass 


Here we briefly comment on the case where the Hubble mass term exists in the potential: 

1/(0) = (2.33) 

2 n 

The equation of motion of 0 is given by 


+ 3H(j) + + A0”-^ = 0 


(2.34) 


The equation of motion of the rescaled held (p is the same as (2.6) except that is now 
replaced with 

2 _ 2(6p — 3pn + n) 2 




+ cp . 


(n-2)2 

In particular, for n = nc, the equation of motion becomes 




-36 + c(n + 2)2^^ (n + 2)2 A 

^+7P7Z -= 0 . 


9(n-2)^ 


9(n - 2)2 Hi 


(2.35) 


(2.36) 


As long as |c| 1, the existence of the Hubble mass term does not qualitatively affect the 

pseudo scaling dynamics studied in the previous subsection. If |c| > 1, the crucial difference 
from the case of |c| <C 1 is that the initial position of 0 cannot be taken freely. In particular, 
for positive c > 1, we have 0* = 0 and there is no dynamics. 

In terms of y, it means that the quadratic term dominates the effective potential. 

If p changes, e.g. due to the inflaton decay, the pseudo scaling behavior terminates soon. 
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Let us consider the case of negative Hubble mass: c < — 1. In this case, 0 initially sits 
at the minimum of the potential (p = (pn where 




-cH^^ 


2 \ ^/{n-2) 


(2.37) 


This is the same order of (pmin during inflation, and hence tp begins oscillation around (prain 
with an oscillation time scale of order one in terms of s = ln(f/fj). Note that the temporal 
minimum of the potential (pn exhibits the same scaling as the scaling solution (p cx . 

Thus (p{t) follows the temporal minimum (pnit) with oscillation around it. In this case, this 
oscillation time scale is close to the Hubble time scale. This is qualitatively different from 
the case of |c| 1 in the previous subsection. 

Let us see the dynamics more closely. In order to evaluate the initial amplitude A(p 
around the temporal minimum (pn, we must specify the origin of the Hubble mass term of 
(p. It arises from the Planck-suppressed interaction as 


where ck and cy are constants of order unity, I denotes the inflaton hied and l^nf is the 
inhaton potential. During a slow-roll inhation, the potential energy of the inhaton dominates, 
hence we have c = cy. After inhation, the inhaton begins a coherent oscillation and the 
kinetic energy and potential energy are equally distributed, roughly speaking. In the case 
of harmonic oscillation of the inhaton, we have c = {ck + cv)/‘2. Therefore, unless ck = cy, 
the temporal minimum (pn just after inhation is diherent from that during inhation. This 
shift of the potential minimum occurs in the time scale of inhaton mass, rather than the 
Hubble scale, so (p cannot track the temporal minimum 13 . Thus the (p oscillation around 


the temporal minimum (pu is induced just after inhation, and the amplitude is of the order 
of A0 


(pH 


#10 


Another term is coupling to the Ricci scalar R\ 


- C 


2V 




CR(p‘^ 

6Mj, 




(2.39) 


If this is a dominant origin of the Hubble mass term, c = cr during inhation and c = cr/A 
after inhation for the harmonic inhaton oscillationj^Hl 

Anyway, the oscillation amplitude around the temporal minimum A(p decreases as the 
universe expands according to 


H{A(pY oc a 


-3 


■H- 


AcP oc 


13 


Even ii Ck = cv 1 (p cannot track the temporal minimum 
During the inflaton oscillation, R oscillates between positive and negative values. 


(2.40) 
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because of the number density conservation. By using this, we obtain the following scaling: 

Ar/) 2(nc—n) 

- ~ i _ 

(pH 



The oscillation amplitude relatively decreases compared with the position of the temporal 
minimum for n > Uf., while the ratio is constant for n = ric- Thus 0 safely relaxes to its 
temporal minimum for n > ric- For n = ric, (p oscillates around (pn and never cross cp = 0 
thereafter if A(p < cpn initially. This fact is important to study the formation of topological 


defects, such as domain walls in the axion model 



3 Cosmological implications 


So far we have considered the dynamics of scalar held with n = Uc having only the non- 
renormalizable potential cp"^ and the Hubble mass term. More realistic models may have 
other terms in the potential. In this section we consider cosmological implications of the 
scalar dynamics with n = ric based on the potential 

1/ = ±-mV + -0"- (3.1) 

2 n 

The presence of the Hubble mass term like cH^cp^ does not change the result much as long 
as |c| <C 1. In this type of the potential, the dynamics of (p is eventually dominated by the 
mass term at H < m, even if (p obeys a pseudo scaling law at hrst. This is easily seen by 
rewriting the equation of motion as 


V?" + 


3p — 


n + 2 


n 


cp'+{p^± 


2 2 

pm 

ip~ 


,nx. - n 


(3.2) 


since /i ~ 0{1), the effect of the mass term becomes important sX H < m. Note that even if 
the potential is dominated by the mass term, the scalar held does not always feel the mass 
term if (p obeys a pseudo scaling law as we explained at the end of Sec. 2.2]^^^ 


Fig. I shows numerical results for the time evolution of (p{t) for n = 6 for p = 2/3 with 
a mass term. In the left (right) panel we have assumed positive (negative) mass squared. 
The red solid (green dashed) line corresponds to m = 10“® (m = 10“^^). Other parameters 
and initial conditions are the same as those in Fig. These results clearly show that the 
pseudo scaling law is maintained even after the mass term dominates the potential and the 
oscillation by the mass term starts at at FT ~ m. 


3.1 Abundance and curvature perturbation 

First consider the abundance and the curvature perturbation of (p. Since (p is light during 
inhation, it obtains a quantum huctuation of S(pi ~ iFinf/(27r) and (p may act as a curvaton- 
like held. However, it is stated in Ref. that the curvature perturbation vanishes for both 

In fact, (j){H = m) ~ (pmiPm/Pi) "C dm, where dm is defined below Eq. (|3.5|). 
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Figure 3: Time evolution of (j){t) for n = 6 and p = 2/3. In the left (right) panel we have 
assumed positive (negative) mass squared. The red solid (green dashed) line corresponds to 
m = 10“® {m = 10“^^). 


n > ric and n = n,, because of the scaling behavior of the curvaton held. If 0 obeys a scaling 
solution (2.11), the held value of 0(f) at some uniform density slice does not depend on the 
initial position (pp. this is easily checked as 


0(f) ~ 




(3.3) 


hence the 0j dependence vanishes. Here fi ~ (A0/ ^) is the time at which the ehective 


mass becomes equal to the Hubble scale (see Eq. (2.23)). Thus curvature perturbation is 
not generated by 0 with a scaling law. 

This is true for n > Uc in which the scaling solution is an attractor, but for n = ric we 
have seen that the actual dynamics is far from the scaling law. For fi < f < f 2 , 0 scales as 


<Pit) 


4/(n-2) ^ 


(3.4) 


The initial 0* dependence remains and hence the perturbation 6(j)i is not always erased. 

For concreteness, let ns suppose that 0 decays into radiation at H = with a constant 
decay rate F,^ after the dynamics is dominated by the quadratic mass term for fi < f < f 2 
and oscillates around the origin. We also assume that the inflaton energy density scales 
as Pi cx until the curvaton enters in the quadratic oscillation regime, and decays into 
radiation before the curvaton decays, i.e., we assume F^ < F/ < m with F/ being the inflaton 
decay rate. In this subsection we mainly consider the case of positive mass squared +m? 
and give only results for the case of negative mass squared —m? at the end. For details on 
the case of negative mass squared, we refer readers to App. 0 
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The energy density of the curvaton at the onset of its oscillation, H = m, is evaluated as 


p^{H = m) ~ = mf ~ , (3.5) 

where we have dehned (j)m = {m ?/It represents the boundary at which the mass 
term begins to dominate the potential. Therefore, the fraction of the 0 energy density at 
H = V^ is given by 


pYh = t^) ^ 1 

Pr{H = T^)- Q\Mp) \(j)J \tJ \m) 


This expression assumes that 0 is subdominant at the decay. Otherwise, it is equal to 1. 

From the 


Now, let us evaluate the curvature perturbation generated in this scenario 


#13 


expression (3.5), we can see that the initial 0* dependence explicitly remains. Thus the 
curvature perturbation of the order of C ~ ^(pi/Yi is generated if 0 eventually dominates the 
universe. 

According to the 6N formalism 


14 , the large scale curvature perturbation ( is given 


by the fluctuation of the e-folding number N{x) between the initial spatially flat surface 
where the observable scale exit the horizon, and the hnal uniform density surface which we 
can take H = The total e-folding number can be conveniently divided into N{x) = 
N^^\x) + N^'^Yx) + N^^\x) with 


= In 


— i/end) 
a{H = \ ’ 


Ar(2) = In 


a{H = Ti) 
a(Tr = hfend)J ’ 


= In 


a{H = T^) 
a{H = Tj) 


(3.7) 


where ifjnf and H^nd are the Hubble parameters at the horizon crossing and the end of 
the inflation, respectively. The curvature perturbation is expanded in terms of the inflaton 
fluctuation 6Ii and 50* as 


C(f) ^ dN{x) XX Nt6Ii{x) + N^54>i{x) + 


(3.8) 


Note that the inflaton mostly contributes to A^i and 0 can have its largest contribution to 
A 3 , hence we can safely approximate as Nj ~ Nj^\ A^ ~ A® and A^^ ~ Thus we 

immediately have 



V)nf 

MUdiViniY 


(3.9) 


To derive and A^^ we use the relation at the curvaton decay surface H{x) = T.^ 


15 


Pr{H = T^,x) - 1 - pYH = T^,x) = ST^Mp, (3.10) 

Here we concentrate on the large scale curvature perturbation. 
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where pr denotes the radiation energy density that is created by the inflaton decay. They 
are expressed by using the e-foldings of the scale factor 

Pr{H = T^,x) = pi{H = 

p^{H = r^,x) = p^{H = 

By differentiating this expression with respect to <pi, we obtain 


( 3 . 11 ) 

( 3 . 12 ) 


^(3) _ ^ ^<piP<l> 


= r 


dl 




where 


r = 


3p0 

^P<P 

Apr -\- ‘Pp<f> 


(^2 _|_ 2^^ 


^ct>i P4> 

^P4> 


(3.13) 


(3.14) 


H=r 




with the over-lines indicating the averaged quantities. 

In the case of positive mass squared +m^, to the second order in Sep, we obtain the 
curvature perturbation as 

2 


2rS(j)i{x) 

C{x) = -7— + r 


2(2r -|- 
9 


S(j)i{x) 


(3.15) 


where we have neglected the inflaton fluctuation. The non-linearity parameter of the local 
type is given by, 


6 , 9 

5^"" - ^ 


2(2r -|- r^) 


9 


m 


(Nf+N^^r 


( 3 . 16 ) 


Here we also show the results in the case of negative mass squared —If the mass 


term dominates the dynamics during the regime t 2 < t < the curvature perturbation is 
given by 


C = Nj6I + {n + l)prf3 


1 f 6(j)i 


‘ <<>T 2 ' ’ 


( 3 . 17 ) 


( 3 . 18 ) 


where P ~ 0(10 ^) is a numerical factor. The non-linearity parameter is given by 

6 1 

5 ^^^ {n + l)pr/3{Nf + Niy' 

For more details, see App. 0 

This is one example of the property of the curvature perturbation from the pseudo scaling 
curvaton. Although these expressions depend on how the background evolves and when the 
curvaton decays, it generally holds that the curvature perturbation is generated for n = ric 
in which (p obeys a pseudo scaling law. 


The small scale perturbations do not grow even in the presence of the tachyonic mass if the momentum 
is larger than the tachyonic mass. 
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3.2 Domain wall problem 


Next, we consider implications on the domain wall formation. Consider the case of negative 
mass squared: —The potential has a symmetry 0 —>■ —0 if n is even. If 0 obtains 
a large initial value during inflation and n < ric, it oscillates around the symmetric point 
0 = 0 after inflation and large fluctuations of 0 particles (and other particles which interact 
with 0) are induced due to the efficient particle production. Such large fluctuations may 
cause the non-thermal restoration of the symmetry and 0 may be trapped at the origin 
(0 = 0) 17 . The trapping effect becomes diluted as the universe expands, and 0 hnally 


rolls down to the minimum of the potential. Eventually, the domain wall is formed and the 
universe is expected to be overdosed. This is the rough picture of the domain wall problem 
we consider here. A physically well-motivated example is the Peccei-Quinn (PQ) held, which 
often leads to the axionic domain wall problem 18 due to the nonthermal restoration of the 


PQ symmetry 19,20 


On the other hand, for n = ric, the scalar held 0 obeys the pseudo scaling law. As 
we saw earlier, the energy density of (p is negative for a reasonable initial condition (see 
Eq. ( 2.18| )). Hen ce 0 never crosses the origin during the pseudo scaling and no particle 


production occurs, Moreover, even after the dynamics is dominated by the mass term. 


0 never crosses the origin because the mass term also prevents p from crossing the origin 
(see Eq. (3.2)). The held 0 just oscillates around 0m and the domain wall problem does not 
occur. Therefore, the pseudo scaling solution provides a natural solution to the domain wall 
problem. 

Here we emphasize that the pseudo scaling solution for n = n,, is more powerful in 
solving the domain wall problem than the scaling solution for n > ric- To see this point, let 
us consider the following potential: 


V 


2 

m 


+ h 

i n 


(3.19) 


where we assume n,nc > I and I > 2. First consider the scaling solution n > ric- In this 
case, if the initial value of 0 is large enough, the held 0 obeys the scaling solution after 
inhation. If 0; = S> 0m; the dynamics is dominated by the 0^ term and 0 starts 

to oscillate around the origin when 0 ~ 0^. Thus, for sufficiently large 0/, the domain wall 
problem occurs even if 0 obeys the scaling solution at hrst. 

On the other hand, in the case of the pseudo scaling solution n = Uc, the dynamics 
becomes dominated by the 0^ term when H Hi where Hi is dehned as 


#16 


I!l 

Vx 


n — 2 

A A 2(n-0 

V 


n — 2 



(3.20) 


Even if > 0, no efficient particle production is expected to occur since the period of the oscillation 
is much longer than the cosmic time in the case of the pseudo scaling solution. 

Note that Hi/m = (0//0m,)^"~^^^^, and hence the dynamics is dominated by the b* term before the 
mass term becomes effective for 0; ^ 0m- 
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We can deduce this as follows. In terms of 99 , the effective potential is given by 

2(n-Z) 


1/2 2 I 


Ves{p) ~ 2 y 7f2 




y + 


nH? 




(3.21) 


The amplitude of y is of the order of (^min, and hence the term becomes effective once the 
condition 


2{n-l) 

X] (HiX ri-2 


m \ H 




0{l), 


(3.22) 


is satished. This condition is equivalent to Eq. (3.20). Therefore, for example, if the <p^ term 


dominates the dynamics for the regime ti < t < t 2 , the held value oi (j) at H = Hi is given 
by 


<P{H = TT,) ~ ( A 


(3.23) 


If the initial amplitude satishes (pi ^ (pi, (j) starts to oscillate with the amplitude much 
smaller than (pi. Thus, the domain wall problem is less likely to occur even in this case 
compared to the scaling solution. 

So far we have focused only on the domain wall problem, but similar implications can be 
obtained for other types of topological defects. 


4 Conclusions and discussion 

We have shown that a scalar held with potential 0" with n = ric exhibits a peculiar behavior 
which we call a pseudo scaling. It resembles neither oscillating nor scaling solution. We have 
also shown that the primordial curvature perturbation is generated by this pseudo-scaling 
curvaton. The critical exponent ric depends on the background expansion law: Uc = Q 
{ric = 10) in the MD (RD) universe. Even a renormalizable term can lead to a pseudo- 
scaling dynamics, i.e., n < 4 is obtained for p > 1, where the Hubble parameter is given by 
H = p/t. 

This observation has implications on some realistic models and cosmological scenarios. 
In supersymmetric (SUSY) theories, for example, it is often the case that a scalar held has a 


non-renormalizable potential. An example is the Afheck-Dine (AD) held 10 which consists 


of the hat direction in the MSSM and have the superpotential of the form 11,21 


W = 


7A 


where (p denotes the AD held and the scalar potential is then given by 


V = -h h.c.^ -f- 7 ^| 0 A^ 


(4.1) 


(4.2) 
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where m and Am represent the effect of SUSY breaking. If there is no Hnbble mass term or 
its coefficient is small, and the last term dominates the potential, the AD held may enter the 
psendo-scaling regime. The LH^ hat direction, which is lifted up by the £ = 4 superpotential, 
matches the condition for the pseudo scaling if the background is matter-dominated as in the 
inhaton oscillation dominated era. In such a case, the lepton or baryon asymmetry generated 
by the AD mechanism will be signihcantly suppressed compared with the naive expectation, 
e.g., the case of negative Hubble mass term. 

Another example is the SUSY axion model, in which the superpotential is given by 


W = 




( 4 . 3 ) 


where 0 and 0 are the PQ held having PQ charge -|-1 and —£ respectively. The scalar 
potential is 


V = — -I- 7^ 


l2£ 


(4.4) 


where we have taken 0 = 0 for simplicity. If > 0, 0 can have a large VEV which may be 
suitable as a solution to the strong CP problem [9,23 


In the case where the PQ symmetry is broken before inhation, the axion held usually 
obtains a large huctuation and may cause the isocurvature perturbation problem. One way 
to avoid the isocurvature problem is to consider the situation where the PQ held takes a 
very large value during inhation 24 . In such a case, the isocurvature perturbation gets 


suppressed due to the large expectation value of the PQ held. However, since the initial held 
value is large, particles are expected to be ehectively produced due to the oscillation of the 
PQ held after inhation. As a result, the PQ symmetry may be non-thermally restored and 
hence the domain wall problem may occur after the QCD phase transition 19,20 . In this 
paper, we have seen that if the PQ scalar obeys the pseudo scaling law after inhation, no 
such a disastrous situation is realized. Thus, for example, if / = 3 and the PQ scalar starts 
to move during the inhaton oscillation dominated era (with the mass term), the domain wall 
problem as well as the isocurvature perturbation problem can be naturally solved. 
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A Expansion law in the inflaton oscillation dominated 
era 

In this appendix we derive the expansion law a oc or H = p/t in the inflaton oscillation 
dominated nniverse. We assume that the inflaton potential around the potential minimum 
is approximated as 

V"inf(/) OC P, (A.5) 

where I denotes the inflaton fleld measured from its potential minimum. The equation of 
motion of the inflaton fleld is given by 

I + 3Hi + ^ = 0. (A.6) 

From this we obtain 

Pi + 3Hi^ = 0, (A.7) 

where pi = P12 + Idnf(-f) is the inflaton energy density. On the other hand, we have 

^{li) = P-qV^i-3Hli. (A.8) 

at 

Taking time average over the inflaton oscillation time period and assuming that the Hubble 
parameter is much smaller than the inflaton mass scale, we obtain the Virial theorem 


(/"> = 


(A.9) 


where (■ ■ ■) indicates time average over the inflaton oscillation period. By using this relation, 
we obtain 


Pi + 


6q 

q + 2 


Hpi 


0 . 


(A. 10) 


Thus the energy density scales as 

Pi OC a . 


(A.ll) 


On the other hand, the Freedman equation 3H^Mp = pi implies pi oc 
obtain the expansion law 

+p 5' + 2 

a oc r, p= ——. 

3q 


Therefore, we 


(A. 12) 


For g = 2 (g = 4), the usual MD (RD)-hke universe p = 2/3 (p = 1/2) is realized. More 
unconventional expansion laws are also obtained, for example, p = 1 for g = 1 and p = 5/3 
for g = 1/2. In terms of g, the critical exponent ric is given by 


Tie = 2(g + 1). 


(A,13) 


18 



B Formal solution of a pseudo scaling scalar field 


In this appendix, we give a formal solution of the equation of motion (2.14). Note that the 
energy density of ip 


Pip 2^ 

is conserved. Thus, the kinetic term of ip is given by 

ip' = ±^2 (p^ - V^ii{ip)). 

By integrating this, we obtain the formal solution of ip as 


Si = 


dip 


for ip' > 0 region and 


s - Si = 


V2p ■\/p^ 

dip 
J ip \J Pip 


(B.l) 


(B.2) 


(B.3) 




(B.4) 


for <p' < 0 region. Here Si and ipi are some initial values of s and <p, respectively. In particular, 
the period of the oscillation of ip is given by 


r‘P+ 

= V2 

J (Z7_ 


dip 




(B.5) 


where ip± are the solutions of the equation = Ves_{ip±) with > ip_*^'^ Note that ip± 
as well as p^p depends on the initial condition. However, it is difficult to solve this equation 
analytically. Moreover, the approximation given in Sec. 2^ is sufficient for our purpose in 
this paper. 


C Curvaton with hilltop potential 

In this section we briefly review the curvature perturbation generated by the curvaton with 
a hilltop potential of the form 

^(0) =+ ^0’"- (C.l) 

As shown in Sec. if the initial position 0* is set 0m "C 0* < 9?min, where 0m = {Tn?/ 

is the minimum position of the potential and ipmm ~ {Hf /0 obeys a pseudo scaling 

Here we assume that < 0. For > 0, we must replace by zero and also double the right hand 
side of the equation when n is even. 
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law until H becomes comparable to m. For concreteness, let us suppose that this breaking 
of pseudo scaling law happens at t 2 < t < The held value at this stage is given by (2.29): 


in+l 




(C.2) 


and this can be much smaller than (pm- The scalar dynamics thereafter is the same as that 


with a scalar potential (C.l) with initial position 0* close to the origin. In other words, the 


pseudo scaling scalar dynamics naturally provides us with a hilltop initial condition. In the 
following, we treat 0*(<C (pm) as if it is the “initial” condition of the dynamics, although the 
true initial position is (pi.^^^ 

As in Sec. 3.1, let us assume T^ < T/ < m. The abundance of the curvaton at H = V^ 
is evaluated as 


p^{H = T^) ^ n - 2 

Pr{H = T^) 6 


Mf 


m 


T/ 

T^ 


1/2 




3p-2 


(C.3) 


This expression assumes that (p is subdominant at the decay. Otherwise, it is equal to 1. 
Here hfos and (pos denote the Hubble parameter and the curvaton held value at the onset of 
its oscillation, respectively. What is nontrivial in the hilltop scenario is that ihos strongly 
depends on the initial condition 0* in the hilltop limit (p^ (pm 25,26 . 

Let us express iLos and (pos in terms of 0*. Before the onset of oscillation, we can neglect 


the second term of ( |C.l ). The formal solution of the equation of motion with a tachyonic 
mass term —m^(p‘^/2 and H = p/t is written as 


(p{t) = 0*2^r(i/ + l){mt) 


V = 


3p- 1 


(C.4) 


where lv{x) is the modihed Bessel function of the hrst kind. In general, the solution is 
a superposition of the modihed Bessel function of the hrst and second kind. However, 
without loss of generality we can consider the region (p > 0 and (p > 0, and hence the Bessel 
function of the hrst kind is sufficient to describe the dynamics. The solution is conveniently 
approximated as 


(p{t) ~ 


1 + 


{mt) 


2 1 


2(3p + 1) 


mt 


2'^r(z/ + l) e 


for mt 1, 
for mt S> 1. 


(C.5) 


First note that the velocity of the scalar held is given b}^^ 


0(t) ~ 


P 


m 


3p + l H 
m(p 


for mt <C 1, 
for mt 3> 1. 


(C.6) 


Here and in what follows we focus on the region p > 0 without loss of generality. 
#19 por \V''/H'^\ ^ 1, we have a solution for a general potential: Hcj) = —pV'j{3p + 1). 
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Since the oscillation epoch is signihcantly delayed in the hilltop limit 0* -C 0m, i-e., m/Hos 3> 
1, the condition of the start of the oscillation is given by 


i^0. 


~ 1 f-)- ifos ^ 


m0o 


(C.7) 


On the other hand, from Eq. (C.5), we obtain 


ln(0os/0* 


— 


(C.8) 


For reasonable choice of 0* 0m, we have 0 ~ (9(10 0. This justihes our assumption 

Hos -C m. Then following relations are derived, 


0 00 


^200 


0 00 


00 * 1 + 00 *’ 002 (1 + 0)3 02 ■ 

Using these equations and noting that = T^) oc we obtain 

3p0 0^^P0 3p0 


(C.9) 


P4> 


p(j> 


(C.IO) 


where we have taken only leading terms in 0. Hence the curvature perturbation (3.8) is 
given by 


C = Ni5I + pr(3 


1 f 6(j)^ 


(C.ll) 


to the leading order in 0. If 0* would be a true initial condition set during inflation, then 
50* ~ Trinf/(27r) and this expression gives a correct estimate for the curvature perturbation 
from the hilltop curvaton. In this case, the non-linearity parameter of /nl = 5/(6pr0) ~ 
(9(10) for r ~ 1 is a generic prediction 26 . 

In our pseudo-scaling curvaton case, the true initial position is 0, and it is 50* that 
obtains a quantum fluctuation of ~ i9inf/(27r). From the relation (C.2), we obtain 


c 


Ni6I + (n + l)pr0 


^ 1 / 50A^ 
4>i 2 0 0* y 


(C.12) 


The non-linearity parameter is given by 


6 

5 


/nl 


1 Nl 

{n + l)pr0 {Nf + iV|)2 ■ 


(C.13) 


Therefore, if the inflaton contribution to the curvature perturbation is negligible, we have 
/nl r\j j3 ^/n (9(1-10) for r ~ 1. Thus in this case the non-linearity parameter may 

become sizable. 
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